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We introduce a novel mechanism for itinerant ferromagnetism, which is based on a simple two- 
band model, and using numerical and analytical methods, we show that the Periodic Anderson 
Model (PAM) contains this mechanism. We propose that the mechanism, which does not assume 
an intra-atomic Hund's coupling, is present in both the iron group and some / electron compounds. 



Even though itinerant ferromagnetism was the first col- 
lective electronic phenomena studied quantum mechani- 
cally, the microscopic mechanisms driving this phase are 
still unknown (l]j^]. In 1928, Heisenberg || formulated 
his spin model to address this issue, but as Bloch ||] 
pointed out, a model of localized spins cannot explain 
the metallic ferromagnetism observed in Fe, Co and Ni. 
After seven decades of intense effort we still do not know 
what is the minimal model of itinerant ferromagnetism 
and, more importantly, the basic mechanism of ordering. 

In 1963, Hubbard J5|] and others introduced the Hub- 
bard model to explain the ferromagnetic (FM) prop- 
erties of the iron group, incorporating the kinetic en- 
ergy in a single nondegenerate band with an intra-atomic 
Coulomb repulsion U . With the exception of Nagaoka's 
H and Lieb's j?J theorems, subsequent theoretical ap- 
proches were not controlled enough to determine whether 
the Hubbard model has a FM phase. The central issue is 
the precise evaluation of the energy for the paramagnetic 
(PM) phase. Because it does not properly incorporate 
the correlations, mean field theory overestimates this en- 
ergy and predicts a large FM region jij. In contrast, 
numerical calculations have narrowed the extent of this 
phase to a small region around the Nagaoka point j(| . 

Going beyond the simple one-band Hubbard model 
is advocated, for example, by Vollhardt et. al. 0. 
They note that the inclusion of additional density-density 
interactions, correlated hoppings, and direct exchange 
terms favors FM ordering. A very simple analysis shows 
that increasing the density of states T>(E) below the 
Fermi energy Ep and placing Ep close to the lower band 
edge increases the FM tendency. One can achieve this 
by including additional hopping terms. The effectiveness 
of a next nearestneighbor hopping t' was studied numer- 
ically by Hlubina et al || for the Hubbard model on a 
square lattice. They found a FM state when the van 
Hove singularity in T>(E) occurred at Ep. However, this 
phase was not robust against very small changes in t' . 

Years ago, Slater || and van Vleck jio) speculated 
that band degenerancy is an essential precondition for 
itinerant ferromagnetism. They suggested that the intra- 
atomic Hund coupling in open shells could be transmitted 
from one atom to another by the conduction electrons. 



However, there are FM metals like Ni where the influence 
of the Hund's coupling is not clear. In Ni, Hund's cou- 
pling is associated with the 3d 8 configuration which has 
low probability as the main configurations are roughly 
40% of 3d 10 and 60% of 3d 9 |o). The relevant question 
for some transition metals is thus whether a model in- 
volving just the two configurations is sufficient, or are 
other orbitals and Hund's exchange necessary to explain 
ferromagnetism. Furthermore, there are / electron itin- 
erant ferromagnets, like CeRli3B2 PH| , whose only local 
magnetic coupling is Kondo like, i.e. antiferromagnetic. 

The novel mechanism we now introduce emerges from 
a two-band model, such as the periodic Anderson model 
(PAM) . The basic ingredients are an uncorrelated disper- 
sive band hybridized with a correlated and narrow band. 
Missing is an explicit intra-atomic Hund's exchange. We 
show the PAM supports our mechanism by interpreting 
the results of quantum Monte Carlo (QMC) simulations 
with an effective model derived from it. 

We will discuss our mechanism in the context of: 
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and t a hoppings are only to nearest-neighbor sites. When 
t a = 0, the Hamiltonian is the standard PAM. For the / 
electron compounds, the a and b orbitals play the role of 
the / and d orbitals, and t a ~ 0. For transiton metals, 
they correspond to the 3d and 4s orbitals. For U = 0, 
the resulting Hamiltonian Hq is easily diagonalized: 
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2t a cos k Xi for a hypercubic lattice in dimension D. 
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The operators which create quasi-particles in the lower 
and upper bands are: 
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FIG. 1. Illustration of the effective model and the FM 
mechanism. A is the hibridization gap and <5 a is the inter- 
val of energy where the electrons are polarized. 

In Fig. 1, we illustrate the (one-dimensional) non- 
interacting bands for the case of interest: e a close to Ep 
and above the bottom of the b band. If \V\ <C we 
can identify two subspaces in each band where the states 
have either predominantly b (ip subspace) or a (<p sub- 
space) character. The size of the crossover region around 
the points where the original unhybridized b and a bands 
crossed is proportional to | V^/ £& | ; that is, it is very small. 
The creation operators for the Wannier orbitals ip r(J and 
4> Ya associated with each subspace are: 
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where N is the number of sites. The subsets K> and 
K< are defined by: K> = {k : \v,k\ > \vu\} and 
K < = {k : > |uk}- In this new basis: 

H = H+ + Ht = rt^Ur'o + E T tr4l^ 
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Because the U term in involves only the a orbitals, 
the matrix elements of H connecting the <p and ip sub- 
spaces are small compared to the characteristic energy 
scales of the problem (the matrix elements of H within 
the subspaces). To see this we express aJ CT as a function 
of (p\ a and ipta by nrs t inverting Eqs. and p] to find: 
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the <p orbitals have predominantly a character, while 
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the ip orbitals have predominantly b character, aJ CT w 
Er'^r-r'C- E 1^1 « then | Wr | « \W T \. Con- 
sequently, the a subspace becomes invariant under the 
application of H. In addition, because \W \ ^> |W r ^ |, 
we can establish a hierarchy of terms where the lowest 
order one corresponds to a simple on-site repulsion: 
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with U = U\W \ 4 and nf a 
terms, containing three and two Wo factors, are much 
smaller and are essentially the same as the intersite in- 
teractions which in the past were added to the Hubbard 
model to enhance the ferromagnetism jy. Adding H^ff 
to Hq we get the effective Hamiltonian: 

H *ff = E ( T r*-r'V4>r'<x " 
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The ip and cp orbitals form uncorrelated and correlated 
non-hybridized bands: H e ff = + H?. For the (p or- 
bitals we obtain an effective one band Hubbard model 
with the peculiar double shell like dispersion relation 
shown by the thick lines in Fig. 1. 

Particularly for t a = 0, has a very large density of 
states in the lower shell of the <j> band |jj which is located 
near e a . From Fig. 1 it is also clear that the electrons first 
doubly occupy the uncorrelated ip band states which are 
below e a . However, when Ep gets close to e a , i.e. the 
system is in the mixed valence regime, the electrons close 
to the Fermi level go into some of the correlated <p states. 
Then, the interaction term H^ff, combined with the dou- 
ble shell band structure of Hq , gives rise to a FM ground 
state (GS): The electrons close to Ep spread to higher 
unoccupied k states and polarize, which causes the spa- 
tial part of their wave function to become antisymmetric, 
eliminating double occupancy in real space and reducing 
the Coulomb repulsion to zero. The cost of polarizing 
is just an increase in the kinetic energy proportional to 
S a ~ hvpSk, where vp is the Fermi velocity and 5k is the 
interval in k space in which the electrons are polarized. 

To determine the stability of this unsaturated FM 
state, we compare its energy with that of the PM (non- 
magnetic) state. If we were to build a nonmagnetic state 
with only the states of the lower cf> shell, we would find 
a restricted delocalization for each electron because of 
the exclusion of the finite set of band states (k-states) 
in the upper shell. To avoid the Coulomb repulsion U 
for double occupying a given site, the electrons need to 
occupy all k-states. This means they have to occupy the 
<f> states in the upper and lower shells. This restricted 
delocalization is a direct consequence of Heisenberg's un- 
certainty principle, and the resulting localization length 
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depends on the wave vectors, where the original b and 
a bands crossed, that define the size (Ak) of each shell. 
The energy cost for occupying the 4> states in the upper 
shell is proportional to the hybridization gap A. There- 
fore if U is the dominant energy scale in the problem and 
A ^> 8 a , the FM state lies lower in energy than the non- 
magnetic state. Under these conditions, the effective FM 
interaction is proportional to the hybridization gap A. 

This mechanism for ferromagnetism on a lattice is 
analogous to intra-atomic Hund's mechanism polarizing 
of electrons in atoms. In atoms, we also have different de- 
generate (the equivalent of S a is zero) shells separated by 
an energy gap. If the valence shell is open, the electrons 
polarize to avoid the short range part of the Coulomb 
repulsion (again reflecting the Pauli exclusion principle). 
The energy of an eventual nonmagnetic state is propor- 
tional either to the magnitude of the Coulomb repulsion 
or to the energy gap between different shells. The in- 
terplay between both energies set the scale of Hund's 
intra-atomic exchange coupling. 

The FM mechanism just described applies to any finite 
dimension. For a chain of 16 sites we calculated, by the 
Lanczos method, the exact GS of H e tf for e a = —tf,, infi- 
nite U (U ^> \t£\), and different values of V as a function 
of electron concentration n — N e /AN. We found that the 
GS is a nonsaturated fcrromagnet between n = 1/4 (one 
electron per site in the PAM) and n = 3/8. In the local 
momentum regime QMC |l]| and DMRG Q calcula- 
tions report a FM phase in a very similar range of n, in 
contrast to the much broader range found by dynamical 
mean field theory calculations JwJ . The largest magneti- 
zation M is obtained when n is such that the lower shell 
of the <j> band is completely filled and the upper one is 
nearly empty; i.e., when Ep ~ e a . 

In Fig. 2a we plot the PAM's energy per site E/N, 
computed by our QMC method fll2| , as a function of 
the total spin per site S/N for chains of varying length 
N and fixed electron density n Over these chain 

lengths the data collapse, with E(S) /N showing the min- 
imum Eqs/N at a non-zero value of S/N that represents 
a very good estimate of the magnetization M/N for the 
PAM in the thermodynamic limit (TL). In Figs. 2b and 
2c, we show AE/N = (E GS - E{S = 0))/N and M/N 
for one and two dimensional systems as a function of 
l/N. In both cases M/N smoothly varies to non-zero 
values in the TL. The non-smooth variation of AE/N 
in two dimensions is a consequence of shell effects still 
present in a finite sized system, but is clearly sugges- 
tive of its likely extrapolation to a non-negative value 
for very large N. In Fig. 2d we plot our QMC results for 
(E(FM)-E(0))/E(0) as a function of e f for a 8 x 8 clus- 
ter and two different electron densities. As we increase 
e a (starting from below the bottom of the 6-band) AE 
decreases and then increases. The most stable FM state 
occurs when Ep « e as expected from our discussion of 
the effective model. Increasing e a even further, AE ap- 
proaches zero indicating that the spin of the GS decreases 



(see Fig. 2e) and the system becomes a paramagnet. 
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FIG. 2. a) Energy of the PAM as a function of total spin 
for different chain lengths, b) Scaling of the energy difference 
between the FM and the PM GS's of the PAM in ID (same 
parameters as in Fig. 2a) and 2D (V = 0.5, e a = —0.3, and 
U = 4). c) Scaling of the GS magnetization, otherwise, the 
same as b). d) Relative energy difference AE/E(0) vs. e a for 
V = 0.5t b , t a = 0, and V = 0.5, U = 4 and n = 0.29 (circles), 
n — 0.31 (squares), e) Magnetization as a function of e a for 
2D, otherwise, the same as d). 

By changing e a we can evolve the system from the lo- 
calized (e a <C Ep) to the itinerant case (e a « Ep). From 
Figs. 2d and 2e we see that the itinerant and the localized 
FM phases are apparently continuously connected. With 
decreasing e a , \AE\ decreases while the zero tempera- 
ture magnetization M increases. The strong reduction of 
| AS | is a result of the very small effective magnetic in- 
teraction in the localized limit (Jrkky is order V A flT8fl). 
Decreasing e a increases the population of the lower <p- 
shell. Since most of the electrons in the lower (j> shell are 
polarized when e a < Ep, this leads to an increase of M. 
These results are consistent with the observed behavior 
of M and the Curie temperature T c in LazCei-^RhsB^ 
as function of a; 19|. 

Even though we cannot do finite temperature calcu- 
lations with our present version of the QMC method, 
we can discuss, at least qualitatively, the predictions 
of our mechanism for finite temperatures. If we move 
e a above Ep, the number of a-electrons decreases to- 
gether with the magnetization M. A new energy scale 
e a — Ep emerges. We propose that this scale is responsi- 
ble for the finite temperature peak in the magnetization 
of Ce(Rhi_ a; Ru a ;)3B2 |16| that suggests an ordered state 
with high entropy. At T = 0, M is small because of 
the reduced number of a-electrons. When the tempera- 
ture is of the order of e Q — Ep, electrons are promoted 
from the doubly occupied b states to the unoccupied a 
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states which have a large entropy (large density of states) . 
These a-electrons polarize because of the energy consid- 
erations discussed above. The source of the large entropy 
is thus associated with charge and not with spin degrees 
of freedom, which explains why a state with larger M 
has a higher entropy. From this analysis we predict that 
the entropy below T c contains a considerable contribution 
from the the charge degrees of freedom. 

We can also connect our mechanism with the hydro- 
static pressure dependence of T c . To do this we calculated 
\AE\/N by the QMC method as function of increasing t^ 
(Fig. 3a). Here we are assuming that the main effect of 
the hydrostatic pressure is to increase and to leave the 
other parameters unchanged. The order of magnitude 
of \AE\/N, which should be proportional to T c , and its 
qualitative behavior in Fig. 3a are in good agreement 
with the experimental results in CeRhjE^ jyj. We see 
from Fig. 3a that for the itinerant FM case, \AE\/~N is 
of the order of 100°K. This scale is much larger than the 
magnitude of the RKKY interaction |i~7| (~ 1°K) which 
is commonly used to explain the origin of the magnetic 
phase when the a electrons are localized. We also find 
that the FM state appears close to quarter filling and 
disappears for n close to 3/8. 
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FIG. 3. a) Energy difference per site between the FM and 
PM states as a function of tb- b) Influence of the hopping t a 
on the FM state. The lattice size is in unit cells. 

We used a nonzero t a to study the stability of the FM 
phase when S a is varied. This study is also important 
because, in contrast to / electron compounds, in the FM 
transition metals both bands are dispersive. In Fig. 3b we 
see that the FM phase is even more stable for t a ~ — 0.1t b 
than for t a — and becomes unstable for t a ~ 0.05i&- 
The reason for this asymmetric behavior is easy to un- 
derstand in terms of the variation of S a : If t a is negative, 
then the effect of t a on the dispersion of the <j> band is op- 
posite to that of the hybridization V. When t a ~ —0.14 
we get, for the given e a and V, the minimum value for 5 a 
and therefore the most stable FM case. When we depart 
from this value of t a , S a increases, \AE\ decreases, and 
the FM state becomes less stable. 

In summary, we introduced a novel mechanism for itin- 
erant ferromagnetism which is present in a simple two 
band model. The picture just presented, combined with 



our previous results , allows a reconciliation of the lo- 
calized and delocalized ferromagnetism pictures painted 
by Heisenberg || and Bloch Q. The hybridization be- 
tween bands plays a crucial role. We have also consid- 
ered the case relevant for the iron group where the dis- 
persion of the lower band is not negligible. The fact 
that the ferromagnetism is even more stable for finite 
values of t a indicates that our mechanism is relevant to 
explain the ferromagnetism of the transition metals, like 
Ni, where a correlated and narrow 3d band is hybridized 
with the 4s band. It suggests that the ferromagnetism 
in the transition metals can originate, at least in part, 
in the interplay between the correlations and the partic- 
ular band structure, and not solely in the intra-atomic 
Hund's exchange |L(|. In addition, our results explain 
several qualitative features observed in the ferromagnet 
CeRh 3 B 2 0,00. Elsewhere we will discuss the rela- 
tion of our mechanism to RE(Coi_ 2: Si K )2 (RE=Ho, Er) 
and the uranium monochalcogenides (US, USe, and UTe) 
which are ferromagnets [ po[ . 
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